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SECTION  I 
INTRODUCTION 

Until  recently,  consideration  of  aeroelastic  divergence  has  essen- 
tially eliminated  the  forward  swept  wing  as  an  aircraft  design  option. 

The  static  aeroelastic  instability  of  divergence  of  lifting  surfaces  is 
well  known.  Bisplinghoff  [1]  presents  the  classical  trend  of  divergence 
speed  as  a function  of  wing  sweep.  In  the  figure  taken  from  his  text 
(Figure  1),  he  shows  the  divergence  speed  for  a conventional  wing  re- 
duces dramatically  with  moderate  forward  sweep,  but  the  divergence  speed 
becomes  very  high  with  moderate  aft  sweep. 

Bending  deformation  affects  the  aeroelastic  behavior  of  swept 
wings.  For  a slender  wing  with  aft  sweep,  bending  produces  a reduction 
in  the  local  angle  of  attack  known  as  wash-out.  Wash-out  unloads  the 
wing  and  virtually  eliminates  the  problem  of  divergence  in  aft  swept 
wings.  However,  for  a slender  wing  with  forward  sweep,  bending  produces 
an  increase  in  the  local  angle  of  attack,  or  wash-in.  Wash-in  increases 
the  aerodynamic  loading  and  total  flexible  lift  curve  slope  of  the  wing, 
and  consequently  reduces  the  aeroelastic  divergence  speed.  An  approach 
to  the  problem  of  increasing  the  divergence  speed  is  to  reduce  the 
bending  deformation,  and  the  wing  wash-in.  For  the  conventional  metallic 
wing  under  a given  aerodynamic  loading,  the  bending  deformation  is 
reduced  by  increasing  the  wing  bending  stiffness  which  normally  requires 
an  increase  in  structural  material  with  an  associated  increase  in  weight. 
For  a conventional  metallic  wing  structure  with  a forward  sweep  greater 
than  15°,  v^e  weight  required  to  provide  adequate  stiffness  for 
sufficiently  high  divergence  speeds  is  prohibitive. 

A different  approach  to  increasing  divergence  speed  is  possible 
when  advanced  composite  materials  are  used  in  the  wing  structure.  If 
one  looks  carefully  at  the  problem  of  divergence,  only  a reduction  in 
wash-in  is  required,  not  necessarily  an  accompanying  reduction  in  bending. 
Advanced  composites  such  as  graphite-epoxy  and  boron-epoxy  have  signifi- 
cantly higher  specific  stiffness  and  specific  strength  characteristics 
than  conventional  aircraft  metals.  Additionally,  these  properties  are 
directional.  The  directional  properties  of  composites  can  be  oriented 
to  alter  the  deformation  under  loading.  By  orienting  the  composites  in 
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advantageous  directions,  wash-in  of  a forward  swept  wing  can  be  reduced, 
and  hence,  the  divergence  speed  can  be  increased.  Therefore,  to  increase 
divergence  speed,  significantly  less  weight  would  be  required  for  a 
composite  structure  than  for  a conventional  metal  structure. 

The  technology  to  design  for  a desired  aeroelastic  response  of  a 
lifting  surface  using  advanced  filamentary  composite  materials  has  been 
named  aeroelastic  tailoring.  References  2 through  11  describe  the  aero- 
elastic tailoring  technology  and  its  applications.  Krone  “2]  applied 
the  aeroelastic  tailoring  procedure  described  in  Reference  3 to  design 
for  elimination  of  aeroelastic  divergence.  He  showed  that  the  weight  of 
executive  transport  and  lightweight  fighter  wings  with  sweeps  from 
35°  aft  to  35°  forward  could  be  significantly  reduced  using 
tailored  composites.  A weight  comparison  of  a metallic  wing  and  a 
tailored  composite  wing  for  a lightweight  fighter  is  presented  in  Figure 
2 taken  from  Reference  2.  The  figure  shows  that  for  increasing  forward 
sweep  the  weight  required  in  aluminum  to  provide  adequate  stiffness 
increases  much  faster  than  the  weight  required  in  tailored  composites. 
Weisshaar  [5,6]  used  laminated  beam  theory  and  aerodynamic  strip  theory 
to  predict  the  static  aeroelastic  divergence  characteristic  of 
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figure  2.  Lightweight  Fighter  Wing  Skin  Weight  Variation  with  Sweep. 

swept  wings.  He  showed  that,  because  of  elastic  coupling  between  the 
bending  and  torsional  deformation  of  the  wing  box,  laminated  composites 
may  be  used  to  preclude  divergenci  for  a large  range  of  forward  sweep 
angles. 

Since  the  forward  swept  wing  has  not  been  considered  a serious 
design  option,  there  is  a scarcity  of  data  on  its  structural  and  aero- 
elastic  characteristics.  The  Flight  Dynamics  Laboratory  recognized  the 
need  for  experimental  data  that  would  illustrate  the  principle  of  aero- 
elastic  tailoring  with  composites  and  its  application  to  divergence  of 
forward  swept  wings.  This  report  describes  the  design,  analysis  and 
testing  of  an  aer.<Jastic  model  which  incorporates  variable  forward  wing 
sweep.  Four  plates  of  the  same  planform,  one  aluminum  and  three  graph- 
ite-epoxy composite  plates  with  different  laminate  orientations,  were 
individually  incorporated  in  the  model  as  the  structural  element.  The 
test  and  analysis  results  illustrate  a simple,  yet  effective,  form  of 
aeroelastic  tailoring. 
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MODEL  DESIGN  AND  FABRICATION 
1.  DESIGN  CONSIDERATIONS 

The  objective  of  the  model  design  was  not  to  replicate  a full-scale 
wing,  but  rather  to  create  a versatile  research  tool  from  which  a maximum 
amount  of  data  could  be  obtained.  The  relatively  simple  wing  design 
facilitated  fabrication  and  computer  modeling,  yet  it  was  of  sufficient 
complexity  to  provide  experience  in  analyzing  and  wind  tunnel  testing  a 
forward  swept  wing  model. 

The  half-span  model  was  designed  to  diverge  at  approximately  the 
middle  of  the  velocity  range  of  the  Air  Force  Institute  of  Technology 
(AFIT)  five  foot  wind  tunnel.  The  maximum  velocity  of  this  tunnel  is 
approximately  300  feet  per  second.  The  boundary  layer  at  the  maximum 
velocity  is  approximately  three  inches  thick. 

The  model  design  evolved  from  the  basic  concept  of  using  a cantile- 
vered plate  as  the  load  carrying  member  with  an  airfoil  shaped  (NACA 
0010)  polyurethane  foam  sleeve  surrounding  the  plate.  The  model  with 
three  of  the  sleeve  sections  removed  from  the  plate  is  shown  in  Figure 
3.  The  plate  concept  was  required  in  order  to  demonstrate  the  effect  of 
tailoring  of  composite  materials  while  allowing  divergence  of  the  wing 
within  the  available  tunnel  velocity  range.  A conventional  two-skin 
wing  box  design  could  not  be  used  because  it  would  have  presented  an 
excessive  stiffness  problem.  The  half-span  model  was  designed  so  the 
plate  could  be  removed  permitting  testing  of  both  aluminum  and  graphite- 
epoxy  materials  while  using  the  same  aerodynamic  sleeve. 

The  initial  sizing  of  the  plate  was  accomplished  by  using  the 
closed  form  solution  for  divergence  dynamic  pressure  given  in  Bispling- 
hoff  (Reference  1)  for  a uniform  slender  swept  wing: 

6.33EI 

qn  = o 5 (1 ) 

a circos£A|sinA| 

The  lift  curve  slope  aQ  was  assumed  to  be  less  than  the  two-dimensional 
value  (2tt)  for  conservatism. 

The  span  at  30°  forward  sweep  was  24  Inches  with  a full -span  aspect 
ratio  of  4 and  a taper  ratio  of  0.4.  The  aspect  ratio  was  selected 
as  representative  of  current  fighter  aircraft  designs.  For  a wing  of 
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this  geometry,  the  thickness  of  an  aluminum  plate  and  of  a 0°t45"  lami- 
nated graphite-epoxy  plate  was  established  for  a dynamic  pressure  at  the 
midrange  of  the  wing  tunnel. 

To  increase  the  capability  of  the  model,  two  features  were  incorpo- 
rated in  the  model  design.  The  leading  edge  sweep  could  be  varied  from 


Figure  3.  Flight  Dynamics  Laboratory  Forward  Swept  Wing  Model. 


0°  to  60°  forward  in  increments  of  15°,  as  shown  in  Figure  4.  At  zero 
leading  edge  sweep,  the  maximum  span  is  33.0  inches.  In  addition,  the 
ability  to  vary  the  model  angle  of  attack  was  included.  Variable  model 
angle  of  attack  was  required  for  the  subcritical  divergence  projection 
methods  used  during  the  wind  tunnel  testing. 


Figure  4.  The  Five  Forward  Sweep  Positions  of  the  Model. 


2.  COMPONENT  DESCRIPTIONS 

The  structural  load  carrying  member  of  the  wing  model  was  one  of 
four  plates  of  identical  planform.  The  plate  position  and  wing  dimen- 
sions in  the  30°  forward  sweep  position  are  shown  in  Figure  5.  The 

leading  and  trailing  edges  of  the  plate  were  on  the  15  and  65% 

chord  lines,  respectively,  and  the  wing  reference  line  was  on  the  40% 

chord  line.  At  -30°  sweep,  the  plate  tips  were  cut  parallel  to 


the  airstream.  The  pivot  axis  was  located  on  the  reference  line  1.5 
inches  inboard  of  the  wing  root  chord. 

The  4 plates,  1 aluminum  and  3 composite,  were  tested  at  each 
of  the  5 sweep  positions,  allowing  for  information  to  be  collected 
on  20  wing  configurations.  The  0.10  inch,  2024-T6  aluminum  plate 
was  used  as  the  baseline  structure. 

The  3 composite  plates  consisted  of  16  plies  of  NARMCO  T300/ 

5208  graphite-epoxy  with  a nominal  thickness  of  5.25  mil.  Each  of  these 
plates  were  cut  from  a larger  plate  with  a symmetric  layup  of  four  0° 
plies  on  the  outside  of  the  laminate  and  two  pairs  of  t45°  plies  on  the 
inside,  or  [04,(-45,+45)2]s,  (Table  1).  As  shown  in  Figure  6,  one  plate 
was  cut  so  the  0°  plies  were  parallel  to  the  reference  line.  The  second 
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and  third  plates  were  cut  so  the  0°  plies,  and  consequently  the  lami- 
nates, were  rotated  7.5°  and  15°  forward  of  the  reference  line,  respec- 
tively. A positive  fiber  orientation  is  forward  of  the  reference  axis. 
After  curing,  the  16  ply  laminate  had  an  average  total  thickness  of 
0.080  inches. 

After  fabrication,  each  plate  was  instrumented  with  strain  gage 
rosettes  located  on  either  side  of  the  plate  on  the  reference  line  four 
inches  outboard  of  the  root.  The  two  opposing  center  gages  on  each 
plate  were  connected  in  a Wheatstone  bridge  to  record  the  average  bending 
strain.  The  torsional  strain  was  obtained  from  an  opposing  set  of  gages 
oriented  45°  from  the  reference  line.  The  remaining  set  of  45°  gages 
were  used  as  spares.  Since  only  differential  voltage  readings  were 
required  for  the  subcritical  divergence  projection  methods,  the  gages 
were  not  calibrated. 

To  achieve  the  minimum  sleeve  stiffness  and  promote  durability,  the 
sleeve  was  sectioned  and  each  section  employed  a bridging  concept. 

Bridges  are  commonly  used  in  beam  type  flutter  models  to  transfer  loads 
from  the  airfoil  to  the  load  carrying  member  over  a minimum  of  beam 
area.  In  the  midspan  of  each  section,  an  aluminum  U-shaped  channel  was 
encased  in  the  foam  on  both  sides  of  the  plates.  The  aluminum  plate 
with  the  bridges  installed,  prior  to  foaming  of  the  sleeve,  is  shown  in 
Figure  7.  After  foaming,  the  sleeve  was  sectioned  to  reduce  the  bending 
and  torsional  stiffness  attributed  to  the  sleeve.  The  crosswise  dowels 
added  lateral  stability  to  each  section  and  transmitted  the  airloads  on 
the  section  onto  the  bridge.  Only  the  bridges  contacted  the  plates. 

Two  aluminum  bolts  one-half  inch  from  the  leading  and  trailing  edges  of 
the  plate  held  each  section  to  the  plate. 

The  variable  sweep  mechanism  (Figure  8)  cantilevered  the  wing  and 
provided  the  variable  sweep  and  variable  angle  of  attack  features  of  the 
model.  The  structural  plate  root  was  clamped  between  two  0.25  inch 
steel  plates.  The  wing  pivoted  about  a 9/16  inch  bolt  located  near  the 
center  of  this  mechanism.  The  nine  holes  located  on  an  arc  centered  at 
the  pivot  were  used  to  align  the  wing  at  the  desired  sweep  angles. 

A fairing  that  had  been  used  for  a previous  test  in  the  AFIT  tunnel 
was  adapted  for  mounting  the  wing  models.  The  sweep  mechanism  was 
housed  inside  the  fairing  which  was  mounted  to  the  wind  tunnel  ceiling. 
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Figure  7.  Aluminum  Plate  with  Bridges. 

Figure  9 shows  a sketch  of  the  wing  installed  in  the  fairing  at  the  30° 
forward  sweep  position.  The  stainless  steel  tube  on  the  sweep  mechanism 
projected  through  the  fairing  and  tunnel  ceiling,  and  was  used  to  control 
the  wing  angle  of  attack  from  outside  the  tunnel.  The  opening  in  the 
fairing  at  the  wing  root  was  minimized  depending  on  the  wing  sweep  by 
various  cover  plates.  One  side  of  the  fairing  was  hinged  to  provide 
access  to  the  sweep  mechanism. 
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In  order  to  avoid  confusion  in  the  following  sections  of  this 
report,  the  following  terms  will  be  used  to  describe  the  wing  configure- 
tions.  The  term  plate  refers  to  one  of  the  four  plates  without  the 
sleeve,  The  term  rotated  refers  to  one  of  the  three  graphite-epoxy 
plates:  the  nonrotated,  7.5°  rotated  and  15°  rotated  plates.  The  term 
model  refers  to  a plate  and  sleeve  configuration  such  as  aluminum  model 
or  7.5  rotated  model.  The  term  wing  is  a generic  reference  applying  to 
both  plates  and  models. 


SECTION  III 
PRE-WIND  TUNNEL  TESTS 


1.  MATERIAL  PROPERTIES 

As  properties  of  composite  materials  can  vary  significantly  from 
batch  to  batch,  specimens  of  the  NARMCO  T300/5208  graphite-epoxy  were 
tested  in  tension  to  failure  to  experimentally  determined  the  elastic 
moduli  and  ultimate  tensile  strains.  Each  specimen  consisted  of  16 
plies  of  the  graphite-epoxy.  For  the  longitudinal  modulus,  four  speci- 
mens with  all  fibers  oriented  in  the  load  direction  were  tested.  For 
the  transverse  modulus,  four  specimens  with  all  fibers  oriented  perpen- 
dicular to  the  load  direction  were  tested.  For  the  shear  modulus,  four 
specimens  with  ±45°  plies  were  loaded  in  tension.  The  average  measured 
values  for  these  moduli  are  given  in  Table  2 with  the  aluminum  material 
properties. 


TABLE  2 

ELASTIC  CONSTANTS 


Constant 

Aluminum 

Graphite-epoxy 

Density  (psi) 

.1 

.059 

Longitudinal  Modulus  E*  (psi) 

10. 5x1 06 

20.8x10® 

Transverse  Modulus  E^  (psi) 

10.5xl06 

1.54x10® 

Poisson's  Ratio  v12 

.3 

.327 

Shear  Modulus  (psi ) 

Ultimate  Strains: 

4.04xl06 

0.80x10® 

Longitudinal  (yin/in) 

11.2X103 

Transverse  (yin/in) 

4. 7x1 03 

Shear  (yin/ in) 

19.0xl03 

2.  SLEEVE  MASS  AND  INERTIA  DATA 

Each  sleeve  section  was  first  weighed  to  the  nearest  0.1  gram  on  a 
balance.  The  center  of  gravity  of  each  section  was  located  by  balancing 
the  section  on  a knife  edge  at  three  angular  orientations.  Lines  were 
etched  in  the  surface  of  the  section  along  the  knife  edge.  The  inter- 
section of  the  lines  was  the  center  of  gravity  of  the  section. 

The  bifilar  pendulum  method  was  used  to  determine  the  rotational 
inertia  of  each  section  about  an  axis  parallel  to  the  wing  leading  edge 
and  through  the  center  of  gravity  of  each  section.  The  measured  mass 
inertia  and  center  of  gravity  location  of  each  section  are  given  in 
Table  3. 


TABLE  3 

SLEEVE  MASS  DATA 


Section 

Center  of  Gravity 
from  leading  edge 

root  A = 0° 

(in) 

x y 

Mass 

(lb) 

Moment 

of 

Inertia 
About  CG 
(lb-in2 ) 

1 

5.91 

1.94 

.2710 

2.993 

2 

5.44 

5.96 

.2403 

2.338 

3 

4.95 

9.94 

.2077 

1.714 

4 

4.51 

13.96 

.1821 

1.258 

5 

4.02 

18.03 

.1557 

.8874 

6 

3.53 

22.00 

.1299 

. 5838 

7 

3.04 

25.85 

.0931 

.3335 

8 

2.28 

29.30 

.0/88 

.1937 

3.  LOAD  DEFLECTION  TESTS 

All  of  the  plates  and  models  were  loaded  in  bending  and  in  torsion, 
to  determine  the  relative  spanwise  bending  and  torsional  stiffness 
distributions.  Front  surfaced  mirrors  were  imbedded  in  modeling  clay 
and  attached  to  balsa  wood  bridges.  Each  bridge  was  supported  on  the 
surface  of  a sleeve  section  by  two  straight  pins  at  the  trailing  edge 
and  one  pin  at  the  leading  edge.  The  mirrors  were  positioned  on  the 
reference  line  at  the  midpoint  of  each  section.  The  bridges  reduced 
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inaccuracies  that  could  be  caused  by  local  surface  distortion.  A light 
source  was  used  to  reflect  crosshairs  off  the  mirrors  and  onto  a grid- 
board. 

By  recording  the  position  of  the  crosshairs  on  the  gridboard  before 
and  after  loading  the  model,  the  difference  in  twist  or  bending  slope 
between  seccions  could  be  calculated.  The  difference  in  slope,  Ao,  is 
given  by 

AG  = Ah/2 l ' (2) 

where  Ah  is  the  distance  the  crosshairs  move  on  the  gridboard,  parallel 
to  the  reference  axis  for  bending  and  perpendicular  to  the  reference 
axis  for  twist.  The  length,  £,  is  the  distance  between  the  gridboard 
and  the  plane  of  the  wing. 

The  average  bending  stiffness,  El,  between  mirrors  is  given  by 

El  = MAy/AG  (3) 

where  M is  the  average  applied  moment  and  Ay  is  the  distance  between 
adjacent  mirrors.  The  average  value  for  torsional  stiffness,  GJ,  is 

GJ  = TAy/AG  (4) 

where  T is  the  applied  torque.  Stiffness  distributions  for  the  aluminum 
plate  and  model  are  compared  to  the  theoretical  distributions  for  bend- 
ing, Figure  10,  and  for  torsion,  Figure  11.  Except  for  the  low  value  of 
bending  stiffness  at  the  tip  due  to  the  stress  concentration  caused  by 
the  load,  the  measured  plate  stiffnesses  compare  favorably  with  the 
theoretical  values. 

As  described  earlier,  the  plates  are  constant  thickness,  therefore 
the  stiffnesses  vary  linearly  with  span.  Dividing  the  measured  stiff- 
nesses by  the  respective  chord  results  in  a constant  value  along  the 
span  and  provides  a means  of  comparison  of  theory  and  experiment  for  the 
overall  plate.  Neglecting  extreme  deviations  from  the  linear  relation- 
ship, the  results  were  averaged  and  are  given  with  the  symbol  key  in 
Figure  1 0 and  1 1 . 

The  average  plate  bending  and  torsional  stiffnesses  are  both  less 
than  the  theoretical  values,  but  the  difference  is  less  than  1 % . The 
contribution  of  the  model  stiffness  due  to  the  sleeve  is  less  than  6%  in 
bending  but  greater  than  20%  in  torsion. 
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Another  loading  test  was  conducted  to  locate  the  approximate  zero 
twist  axis  of  each  plate.  A perpendicular  point  load  was  applied 
directly  to  the  plate  at  successive  locations  between  two  adjacent 
sleeve  sections  and  the  motion  of  the  reflection  of  the  mirror  on  the 
inboard  section  was  observed.  The  point  where  the  reflection  remained 
stationary  is  or,  the  zero  twist  axis.  Figure  12  shows  the  zero  twist 
axis  for  each  of  the  four  models.  The  aluminum  model  has  a zero  twist 
axis  nearly  perpendicular  to  the  cantilevered  root  line.  The  zero  twist 
axis  is  approximately  oriented  with  the  0°  fibers  on  each  of  the  compos- 
ite model s. 


Figure  12.  Zero  Twist  Axes. 


4.  GROUND  VIBRATION  TESTS  j 

Ground  vibration  tests  were  the  last  tests  conducted  before  entry  { 

1 

into  the  wind  tunnel.  The  frequencies  and  mode  shapes  for  each  plate  | 

and  for  each  model  were  measured.  A roving  accelerometer  was  used  to  | 

record  the  relative  displacement  amplitudes  at  17  points  on  the  wing: 


two  chordwise  locations  on  the  mid-span  of  each  section  at  the  leading 
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TABLE  4 

WIND  TUNNEL  MODEL  EXPERIMENTAL  FREQUENCIES 


Mode 

Aluminum 
Plate  Model 

Nonrotated 
Plate  Model 

Graphi te-epoxy 

7.5°  Rotated  ' 
Plate  Model 

15°  Rotated 
Plate  Model 

IB 

4.40 

3.32 

6.13 

3.57 

5.60 

3.40 

5.18 

3.19 

2B 

21.81 

16.97 

31.29 

17.04 

29.21 

16.96 

23.35 

15.36 

IT 

59.10 

36.09 

41 .53 

31.57 

42.48 

32.04 

47.25 

31  .42 

3B 

60.72 

43.48 

82.05 

46.53 

76.63 

44.36 

67.52 

40.70 

2T 

141.64 

81 .85 

98.60 

64.88 

110.01 

65.88 

113.76 

68.74 

4B 

115.34 

87.07 

152.93 

86.22 

145.77 

83.08 

— 

75.99 

Units:  Hz 
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SECTION  IV 


MODEL  ANALYSIS 

1.  METHODS  OF  ANALYSIS 

Slender  beam  theory  [I]  and  the  method  developed  by  Weisshaar  [5], 
were  used  to  perforin  preliminary  analyses  of  the  four  models  to  insure 
that  divergence  could  be  obtained  within  the  speed  range  of  the  wind 
tunnel.  The  analyses  were  refined  using  two  procedures,  TSO  and  NASTRAN. 
The  analysis  methods  are  described  in  this  section. 

a.  Slender  Beam  Theory 

The  slender  beam  theory  as  presented  by  Bisplinghoff  [1]  assumes  a 
high  aspect  ratio  wing  structure  can  be  modeled  by  uncoupled  bending 
and  torsion  flexibility  coefficients  at  given  span  stations  along  an 
elastic  axis.  The  wing  model  had  a sufficiently  high  span  to  chord 
ratio  and  could  be  analyzed  by  this  method,  but  the  graphite-epoxy 
models  exhibit  coupled  bending  and  torsion  flexibility,  so  only  the 
baseline  aluminum  wing  could  be  analyzed.  Modified  strip  theory  aero- 
dynamics provided  the  aerodynamic  influence  coefficients  required  to 
perform  the  divergence  analyses. 

b.  CUING 

The  computer  procedure  CWING  is  an  inexpensive  analytical  tool 
developed  by  Dr.  T.  A.  Weisshaar  of  Virginia  Tech  for  the  Flight  Dynamics 
Laboratory.  Initially,  the  program  provided  a closed  form  solution  of 
the  divergence  problem  for  specific  wing  configurations  [5].  In  subse- 
quent studies  [6],  Dr.  Weisshaar  expanded  CWING  to  provide  analysis  of 
divergence,  lift  redistribution,  and  aileron  effectiveness  of  more 
general  wing  configurations.  CWING  uses  a Weissenger-L  aerodynamic 
theory  and  a structural  model  that  has  its  properties  defined  at  a 
finite  number  of  span  stations. 

c.  TSO 

The  aeroelastic  tailoring  computer  procedure  TSO  [3,4]  is  an  inter- 
disciplinary preliminary  design  program  combining  aerodynamic,  static 
aeroelastic,  flutter,  and  structural  analyses.  Low  to  moderate  aspect 
ratio  wings  can  be  modeled  as  plates,  therefore  the  direct  Rayleigh- 
Ritz  energy  formulation  for  a plate  is  used  to  perform  structural 
analyses. 
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In  TSQ,  a symmetric  wing  structural  box  is  simulated  by  a trapezoi- 
dal plate  with  depth  and  skin  thickness  given  by  biquadratic  polynomials. 

Three  skin  orientations  may  be  modeled.  In  most  wing  structures,  the 

structural  box  has  a depth  much  greater  than  the  skin  thickness,  and  ), 

hence,  the  TSO  stiffness  polynomial  was  formulated  assuming  the  plies  of 

each  orientation  are  distributed  evenly  over  the  skin  thickness.  Because 

the  wind  tunnel  model  structural  box  was  a plate  with  no  core  between 

upper  and  lower  skins,  the  stacking  sequence  plays  an  important  role  in  i 

the  overall  stiffness  of  the  plate.  To  account  for  the  error  induced  by  * 

the  stacking  sequence  distribution  assumed  in  TSO,  equivalent  thicknesses 

were  calculated  for  a distributed  stacking  sequence  that  would  yield  the  i 

same  flexibility  as  a specifically  distributed  laminate.  The  equivalent  i 

thicknesses  for  each  orientation  were  calculated  by  equating  the  defini-  I 

tion  of  the  flexibility  matrix  for  a distributed  stacking  sequence  given  f 

in  Reference  3 with  the  flexibility  matrix  for  a specifically  distributed 

laminate.  I 

Two  aerodynamic  matrices  used  in  TSO  are  calculated  in  other  comput-  | 

er  procedures.  The  steady  aerodynamic  matrix  is  provided  by  a Woodward  t 

aerodynamic  routine,  ROT  [4,12],  This  matrix  is  used  with  the  structural  j 

influence  matrix  generated  in  TSO  to  calculate  the  divergence  velocity. 

The  unsteady  aerodynamic  matrices  are  provided  by  a doublet  lattice 
aerodynamic  routine,  N5KA  [4,13].  An  aerodynamic  matrix  is  calculated 
for  each  of  20  reduced  frequencies  and  a K-method  modal  flutter 
solution  is  used  to  solve  for  the  velocities,  frequencies,  and  dampings. 

For  a reduced  frequency  near  zero  the  corresponding  aerodynamic  matrix 
approximates  the  steady  aerodynamic  matrix,  and  hence,  the  divergence 
velocities  can  be  calculated  by  the  dynamic  analysis  in  TSO. 

d.  NASTRAN 

The  NASTRAN  finite  element  structural  analysis  computer  program 
[14]  was  used  for  stress  analysis,  free  vibration  analysis,  and  flutter 
and  divergence  analyses.  Levels  16  and  17  were  used  In  all  the  NASTRAN 
analyses. 

For  the  aluminum  plate,  the  homogeneous  elements  CQUAD2  and  CTRIA2 
were  used  with' the  material  properties  input  on  a MAT1  card.  For  the 
graphite  epoxy  plates,  CQUAD1  and  CTRIA1  elements  were  used  to  simulate 
the  anisotropic  properties  of  composite  laminates.  For  each  of  the 
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three  composite  laminates,  representative  3X3  in-plane  and  bending 
stiffness  matrices  were  computed  and  input  on  MAT2  cards,  These  stiff- 
ness matrices  were  obtained  from  program  SQ5  [15],  which  gives  the  in- 
plane and  bending  stiffness  of  a laminate  accounting  for  ply  material 
properties,  ply  thickness,  stacking  sequence  and  orientation. 

A stress  analysis  of  the  three  composite  laminates  was  accomplished 
using  rigid  Format  1 of  NASTRAN.  The  highest  expected  steady  airload 
was  first  obtained  from  computer  Procedure  TSO  using  Woodward  aerodynam- 
ics. These  airloads  were  then  resolved  into  lift  forces  and  moments  at 
the  NASTRAN  grid  points  by  a program  known  as  BEAMING  [16].  The  static 
loads  were  then  applied  to  the  finite  element  NASTRAN  model  to  calculate 
stresses,  element  forces  and  displacements  at  the  grid  points.  The  area 
of  the  inode!  near  the  root  of  the  wing,  where  the  highest  stresses  were 
expected,  was  divided  into  a finer  mesh  for  better  stress  definition. 
Using  the  calculated  element  forces,  point  stress  analysis  program  SQ5 
was  used  to  compute  strain  margins  for  each  ply  of  graphite  epoxy  in  the 
unrotated  laminate  near  the  cantilevered  root  area. 

Rigid  Format  3 of  NASTRAN  was  used  to  extract  the  first  six  normal 
mode  shapes  and  corresponding  -frequencies  using  the  inverse  power  eigen- 
value extraction  method. 

Rigid  Format  10  of  NASTRAN  was  used  for  flutter  analysis.  This 
rigid  format  incorporated  doublet  lattice  aerodynamic  theory  to  compute 
the  aerodynamic  influence  coefficient  matrix  used  in  the  flutter  equa- 
tion. The  K-method  of  modal  flutter  solution  was  used  to  solve  the 
flutter  equation  for  both  flutter  and  divergence  speeds.  A range  of 
reduced  frequencies  down  to  zero  was  used  to  obtain  corresponding  values 
of  damping  and  frequency  at  each  value  of  velocity  for  each  mode  in  the 
analysis.  The  first  three  normal  modes  were  used  in  the  modal  solution. 
Flutter  of  the  wing  was  indicated  when  the  damping  of  the  mode  was  equal 
to  zero.  Divergence  of  the  wing  was  indicated  when  the  damping  and 
frequency  of  a mode  simultaneously  went  to  zero. 

2.  ANALYTICAL  MODELS 

a.  Slender  Beam  Theory  Analysis 

The  beam  theory  divergence  analysis  was  performed  at  the  five  sweep 
positions  for  the  wing  assuming  aerodynamics  for  the  wing  planform  and 


stiffness  due  to  only  the  aluminum  plate.  The  wing  planform  was  divided 
into  seven  strips  of  equal  width.  Theoretical  values  of  bending  and 
torsional  stiffness  were  calculated  at  the  midpoint  of  each  of  the  seven 
strips.  The  aerodynamic  strip  theory  was  modified  by  calculating  the 
aerodynamic  center  and  local  lift  coefficient  for  each  of  the  seven 
strips  using  doublet  lattice  aerodynamic  theory.  The  calculated  diver- 
gence dynamic  pressures  confirmed  that  the  model  would  diverge  well 
within  the  range  of  the  tunnel. 

The  effect  of  the  increase  in  stiffness  due  to  the  sleeve  is  evident 
in  Table  5.  The  divergence  dynamic  pressures  presented  in  this  table 
were  recalculated  using  the  measured  bending  and  torsional  stiffnesses 
of  the  aluminum  plate  alone  and  of  the  aluminum  model.  Theoretical 
aerodynamics  for  the  wing  planform  were  used  in  the  calculation  of  these 
dynamic  pressures. 

TABLE  5 

BEAM  THEORY  DIVERGENCE  DYNAMIC  PRESSURE  PREDICTIONS 
FOR  THE  ALUMINUM  MODEL 


Sweep 

PI  ate 

Model 

0 

28.96 

34.75 

-15° 

14.02 

15.53 

-30° 

11.46 

12.39 

o 

LO 

1 

12.27 

13.10 

-60° 

17.52 

18.57 

Model  aerodynamic  planform  used. 

Units:  psf 


b.  CWING  Analysis 

The  early  version  of  CWING  was  modified  for  use  as  a subprogram  of 
an  analytical  procedure  developed  for  this  effort.  The  main  program 
varied  ply  orientation,  stacking  sequence,  and  sweep,  and  called  CWING 
to  calculate  the  divergence  dynamic  pressures.  This  analysis  was  quali- 
tative since  CWING  analyzed  wing  structures  with  similar  cross  section 
along  the  span,  but  the  plates  had  constant  thicknesses.  For  this  anal- 


V -M  ysis,  wing  stiffness  was  due  to  the  plate  alone  and  aerodynamics  were 


calculated  for  the  entire  planform.  The  divergence  dynamic  pressures 
for  all  of  the  stacking  sequences  analyzed  are  found  in  Appendix  B. 

The  CWING  analysis  showed  that  placing  the  0°  plies  farthest  from 
the  wing  centerplane  resulted  in  the  greatest  divergence  dynamic  pres- 
sures. Thus,  the  [0^, (-45,+45)2ls  stacking  sequence  was  chosen  as  a 
baseline  laminate.  Table  6 presents  the  divergence  dynamic  pressure 
results  of  the  CWING  analysis  for  variations  of  this  lamipate.  The  top 
half  of  Table  6 illustrates  the  increase  in  divergence  dynamic  pressures 
due  to  rotating  only  the  0°  plies. 

TABLE  6 

CWING  DIVERGENCE  DYNAMIC  PRESSURES  FOR  VARIATIONS 
OF  THE  [04,(-45,+45)2]s  GRAPHITE-EPOXY  LAMINATE 


Ply  Angles 
(degs) 

0° 

-15° 

Sweep 

-30° 

-45° 

-60° 

-5, +45, -45 

4.39 

4.03 

4.54 

6.44 

13.66 

0,+45,-45 

8.27 

6.64 

6.82 

9.06 

18.36 

5, +45, -45 

45.43 

14.67 

11.17 

12.59 

23.49 

10, +45, -45 

****** 

314.94 

19.67 

15.93 

26.61 

15, +45, -45 

****** 

****** 

32.80 

16.69 

25.56 

20, +45, -45 

****** 

****** 

36.17 

14.35 

21.49 

25, +45, -45 

****** 

****** 

24.46 

11.03 

16.94 

-5, +40, -50 

4. .53 

4.14 

4.64 

6.57 

13.90 

0,+45,-45 

8.27 

6.64 

6.82 

9.06 

18.36 

5, +50, -40 

35.03 

13.66 

10.79 

12.36 

23.22 

10, +55, -35 

****** 

91.40 

18.24 

15.64 

26.51 

15, +60, -30 

★*★*★★ 

★★*★*★ 

30.59 

16.97 

26.27 

20, +65, -25 

****** 

****** 

39.04 

15.39 

22.95 

25, +70, -20 

****** 

****** 

30.86 

12.37 

18.64 

Units: 

psf 

****** • 

qD  > 103  psf 

The  lower  half 

of  Table 

6 demonstrates  an  increase  in 

divergence 

dynamic  pressures  if 

the  whole  laminate  is  rotated. 

The  7. 

5'“  and  15° 

rotated  laminates  were  chosen 

since  cutting  all  the 

composite  plates  out 

of  one  large  plate  was  more  economical 

than  constructing  plates  of  three 

different  laminates. 

In  addition,  the 

variation  of 

material  properties 

was  minimized.  The 

stacking 

sequence 

and  rotations 

were  described  in 

Table  1 and  Figure  8. 
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c.  TSO  Analysis 

A sketch  of  the  TSO  analytical  model  is  shown  in  Figure  13.  The 
trapezoidal  structural  plate  in  TSO  requires  parallel  root  and  tip 
chords  so  the  tip  was  modified  as  shown.  The  dashed  lines  in  the  sketch 
represent  the  airfoil  planform. 


Results  from  the  material  properties  test  and  a dimensional  check 
of  the  finished  plate  provided  the  required  input  data  of  the  plates. 

The  load  deflection  tests  performed  on  the  plates  were  repeated  in  the 
TSO  analysis  demonstrating  a fairly  accurate  representation  of  the 
structural  plates. 

A comparison  of  the  first  six  modal  frequencies  measured  experimen- 
tally with  those  calculated  analytically  showed  the  analytical  values 
were  higher.  Since  the  analytical  model  was  cantilevered  at  the  root, 
the  lower  experimental  frequencies  were  attributed  to  root  flexibility. 
To  model  the  root  flexibility,  the  analytical  model  was  altered  slightly 
by  moving  the  root  of  the  plate  inboard  (one  inch  for  the  aluminum  plate 
and  one  and  one  half  inches  for  the  graphite  epoxy  plates).  The  analyt- 
ical frequencies  for  the  plates  are  shown  in  Table  7. 

In  order  to  account  for  the  mass  of  the  airfoil  sections  in  the  TSO 
analytical  model,  each  airfoil  section  mass  was  divided  into  three 
lumped  masses.  The  section  masses  and  locations  for  the  TSO  analysis 
are  given  in  Table  8.  The  locations  of  the  masses  are  shown  relative  to 
the  plate  in  Figure  13. 


TABLE  7 


TSO  ANALYSIS  PLATE  NATURAL  FREQUENCIES 


Mode 

A1 uminum 

Graphite-epoxy 

Nonrotated  7.5°  Rotated 

15°  Rotated 

IB 

4.44 

5.98 

4.91 

5.64 

2B 

23.06 

31 .13 

25.27 

29.13 

IT 

59.11 

38.95 

46.66 

41 .94 

3B 

62.70 

84.37 

68.99 

78.99 

2T 

143.90 

98.93 

116.04 

106.17 

4B 

126.26 

168.23 

142.20 

158.90 

Units:  Hz 

TABLE 

8 

TSO  ANALYSIS  SLEEVE 

MASSES 

AND  LOCATIONS 

X 

y 

Wei ght 

X 

y 

Weight 

(in) 

(in) 

(lb) 

(in) 

(in) 

(lb) 

2.60 

1.94 

.1221 

2.39 

5.96 

.1115 

2.17 

9.94 

.0983 

1.95 

13.96 

.0863 

1.73 

18.03 

.0760 

1.52 

22.00 

.0647 

1.31 

25.85 

.0493 

.28 

29.30 

.0281 

6.51 

1.94 

.0685 

5.97 

5.96 

.0527 

5.43 

9.94 

.0416 

4.88 

13.96 

.0326 

4.33 

18.03 

.0221 

3.79 

22.00 

.0155 

3.27 

25.85 

.0053 

1.96 

29.30 

.0075 

10.42 

1.94 

.0804 

9.55 

5.96 

.0762 

8.68 

9.94 

.0678 

7.81 

13.95 

.0632 

6.93 

18.03 

.0576 

6.07 

22.00 

.0497 

5.23 

25.85 

.0385 

3.65 

29.30 

.0432 

The  beam  element  feature  of  TSO  was  used  to  account  for  the  stiff- 
ness added  to  the  model  by  the  sleeve.  This  feature  allows  bending  and 
torsional  rigidity  constants  to  be  input  to  model  linear  spars  and  ribs. 
To  simplify  assigning  values  of  rigidity,  it  was  assumed  that  the  sleeve 
sections  could  be  modeled  by  a pair  of  crossed  beam  elements  with  bending 
stiffness  and  no  torsional  stiffness.  The  locations  are  shown  in  Figure 
13.  Through  iterations  of  values  of  bending  stiffness,  the  first  three 
natural  frequencies  of  the  analytical  model  were  matched  to  the  frequen- 
cies of  the  wind  tunnel  model.  The  locations  and  stiffness  values  of 
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the  Seai!!  are  given  in  Table  9.  T he  final  analytical  natural  frequencies 
are  given  in  Table  10.  The  first  three  natural  mode  shapes  of  each 
model  are  shown  in  Appendix  A. 

TABLE  9 


TSO  ANALYSIS  SLEEVE  BEAM  ELEMENTS  AND  LOCATIONS 


Leading-edge 

Trail ing-edye 

Bending 

Stiffness 

Endpoint 

Endpoint 

El  x IQ*3 

(lb-in2) 

x y 

x y 

(in)  (in) 

(in)  (in) 

A1 uminum 

Graphite 

0.  1.75 

12.46  2.25 

131.75 

240.00 

0.  2.25 

12.46  1.75 

131.75 

240.00 

0.  5.75 

11.38  6.25 

105.40 

192.00 

0.  6.25 

11.38  5.75 

105.40 

192.00 

0.  9.75 

10.30  10.25 

79.05 

144.00 

0.  10.25 

10.30  9.75 

79.05 

144.00 

0.  13.75 

9.21  14.25 

52.75 

96.00 

0.  14.25 

9.21  13.75 

52.75 

96.00 

0.  17.75 

8.13  18.25 

26.35 

48.00 

0.  18.25 

8.13  17.75 

26.35 

48.00 

0.  21.75 

7.04  22.25 

10.54 

19.20 

0.  22.25 

7.04  21.75 

10.54 

19.20 

0.  23.75 

5.96  26.25 

5.27 

9.60 

0.  26.25 

5.96  25.75 

5.27 

9.60 

0.  29.75 

5.43  30.25 

1.054 

1.92 

j 

• 

0.  30.25 

5.43  29.75 

1.054 

1.92 

i 

l 

■ - » 

j i 

I 1 

TABLE  10 

I'/.i 

i " i 

i , ' : 

TSO  ANALYSIS 

WIND  TUNNEL  MODEL 

NATURAL  FREQUENCIES 

i ••  : 
l ■'*>. ; 
i 

Graphite-epoxy 

i - .■ 

( "V 

!,■  • 

Mode 

A1 uminum 

Nonrotated 

7.5°  Rotated  15°  Rotated 

* \-‘  * 

V 1 

IB 

3.30 

3.52 

3.38 

3.11 

r:  i 

| J;  i 

2B 

16.82 

17.78 

16.97 

15.44 

ivj* 

i [ H'. 

IT 

36.15 

31 .36 

31.65 

32.81 

■ ';'*y 

3B 

45.06 

46.75 

42.47 

40.32 

■ *0. 

2T 

83.70 

67.14 

49.11 

61 .51 

4B 

90.18 

83.91 

79.85 

78.63 

■ -1 

Units:  Hz 

1 1 
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After  the  pretunnel  test  data  was  incorporated  in  the  analytical 
model,  divergence  and  flutter  calculations  were  performed.  The  aerody- 
namic paneling  for  the  Woodward  and  doublet  lattice  30°  forward  sweep 
analyses  are  shown  in  Figure  14  and  15.  Similar  paneling  was  developed 
for  the  four  other  sweeps.  The  divergence  dynamic  pressure  predictions 
are  presented  in  Table  11  for  both  the  static  and  dynamic  aeroelastic 
calculations.  The  flutter  dynamic  pressures  and  frequencies  calculated 
in  the  dynamic  analysis  are  given  in  Table  12.  The  flutter  mode  is  a 
coupling  of  second  bending  and  first  torsion.  Frequency  and  damping 
versus  velocity  curves  are  presented  in  Appendix  C. 

TABLE  11 

TSO  ANALYSIS  DIVERGENCE  DYNAMIC  PRESSURES 


Graphite-epoxy 

S pPn  Aluminum  Nonrotated  7.5°  Rotated  i b ' l<o ui  v-h! 

w ^ Static  Dynamic  Static  Dynamic  Static  Dynamic  Static  0yno,.ii 


0° 

42.9 

37.6 

49.7 

40.6 

— 

192.6 

— 

CO 

.V 

-15° 

16.8 

15.9 

15.2 

14.2 

19.2 

17.6 

21  .7 

20.1 

-30° 

13.0 

12.6 

11.3 

11 .0 

12.2 

11.9 

11.7 

11.3 

-45° 

13.1 

12.1 

11.2 

10.4 

11.5 

10.5 

10.2 

9.5 

O 

O 
i O 

1 

16.6 

13.9 

14.1 

11.8 

13.8 

11.6 

11 .9 

9.9 

Units : 

: p-'.f 

A 

comparison  of  the 

divergence  speeds  from 

the  static 

aero  el 

astic 

analysis  with  those  from  the  dynamic  aeroelastic  analysis  (Table  11) 
indicates  some  significant  differences  in  the  divergence  predictions. 
Both  analyses  use  identical  structural  models,  therefore  the  difference 
must  lie  in  either  the  aeroelastic  eigenvalue  equations  or  in  the  aero- 
dynamic analyses. 

A steady  aerodynamic  Influence  coefficient  matrix  can  be  calculated 
in  a doublet  lattice  analysis  by  choosing  a reduced  frequency  near  zero. 
This  steady  aerodynamic  matrix  was  substituted  for  the  Woodward  aerody- 
namic matrix  in  the  TSO  analysis  in  order  to  determine  how  much  differ- 
ence exists  between  the  static  and  dynamic  calculations  of  divergence 
dynamic  pressures.  Five  wing  configurations  were  reanalyzed,  The 


27 


ROOT  FAIRING 


Panel  Divisions 

Chordal se 

Spanwise 

0. 

-.500 

.0625 

-.375 

.125 

-.250 

.1875 

-.125 

.250 

0. 

.375 

.125 

.500 

.250 

.625 

.375 

.750 

.500 

.875 

.625 

1.000 

.750 

.875 

1.0 

Figure  14.  Woodward  Paneling. 
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Chordwlse 

Spanwise  | 

0. 

Ml 

.25 

.180 

.5 

.341 

.75 

.485 

1.00 

■B 

Figure  15.  Doublet  Lattice  Paneling. 


comparisons  between  the  static  and  dynamic  divergence  analyses,  given  in 
Table  13,  show  very  little  difference  in  the  two  eigenvalue  solutions. 
Therefore,  the  difference  in  divergence  results  is  due  to  the  aerodynamic 
methods  used, 

TABLE  13 

COMPARISON  OF  TSO  STATIC  AND  DYNAMIC 
DIVERGENCE  CALCULATIONS 


Wing 

Divergence  Dynamic 
Pressure  (psf) 

Plate 

Sweep 

Static 

Dynamic 

Aluminum 

0" 

37.7 

37.6 

Graphite-epoxy 
7.5°  Rotated 

-15° 

17.8 

17.6 

Graphi te-epoxy 
15°  Rotated 

-30° 

11.4 

11.3 

Graphite-epoxy 
Non rota  ted 

-45° 

10.4 

10.4 

Graphite-epoxy 
7.5°  Rotated 

1 

CT» 

O 

o 

11.8 

11.6 

Two  attempts  were  made  to  resolve  the  differences  in  the  aerodynam- 
ics. The  first  was  in  the  location  of  the  panel  center  of  pressure. 

The  Woodward  analysis  assumes  a center  of  pressure  located  at  the  panel 
midchord,  and  the  doublet  lattice  analysis  assumes  a quarter  chord 
center  of  pressure  location.  For  a fine  paneling  grid,  the  location  of 
the  center  of  pressure  should  have  no  effect.  However,  the  grids  used 
in  these  analyses  are  fairly  coarse.  The  aerodynamic  analyses  were 
performed  at  a Mach  number  of  0.13.  For  speeds  in  this  range,  the  flow 
is  such  that  the  center  of  pressure  should  be  near  the  quarter  chord, 
therefore  the  Woodward  routine  was  altered  to  perform  calculations  with 
the  center  of  pressure  at  the  quarter  chord. 

By  moving  the  center  of  pressure  forward  the  twisting  moment  due  to 
the  same  load  will  be  greater  and,  hence,  the  divergence  speeds  will  be 
lower.  The  TSO  analysis  was  performed  once  again  with  the  recalculated 
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steady  aerodynamic  matrix  on  the  same  wing  configurations  used  to  compare 
the  eigenvalue  problems.  Table  14  has  a comparison  of  the  divergence 
speeds  for  each  of  the  five  cases.  As  expected,  the  divergence  speeds 
calculated  by  the  static  aeroelastic  method  did  decrease. 

The  second  attempt  at  resolving  the  difference  in  the  aerodynamics 
was  to  increase  the  fineness  of  the  Woodward  paneling  without  altering 
the  center  of  pressure  location.  The  number  of  panels  shownJn  Figure 
14  was  doubled,  the  aerodynamics  recalculated,  and  divergence  analyses 
performed.  The  results.  Table  14,  show  a decrease  in  the  divergence 
speeds  comparable  to  the  results  obtained  by  changing  the  center  of 
pressure  location. 

TABLE  14 

CHANGE  IN  TSO  DIVERGENCE  CALCULATIONS  DUE 
TO  WOODWARD  CENTER  OF  PRESSURE  LOCATION 
AND  AERODYNAMIC  PANELING 


Wing 

Plate  Sweep 


Aluminum 


7.5°  Rotated 


15°  Rotated 


Nonrotated 

Graphite-epoxy 
7.5°  Rotated 


Divergence  Dynamic  Pressure  (psf) 

Cp  at  .5c  Cp  at  ,25c  160  Panels  Dynamic 


0° 

42.9 

40.7 

40.5 

37.6 

-15° 

19.2 

18.6 

18.6 

17.6 

-30° 

11.7 

11.5 

11.4 

11.3 

-45° 

!1 .2 

10.8 

11.0 

10.4 

O 

O 

to 

t 

13.8 

13.1 

13.5 

11.6 

d.  NASTRAN  Analysis 

The  NASTRAN  model  used  for  stress  analysis  is  presented  in  Figure 
16.  The  highest  steady  airload  expected  during  testing  was  applied  to 
the  model.  This  airload  condition  occurs  at  0°  sweep,  3°  angle  of 
attack,  and  80%  of  the  predicted  divergence  velocity  for  the  unrotated 
composite  laminate.  A maximum  strain  of  one  seventh  of  the  ultimate 
strain  was  predicted  in  the  plies  located  near  the  root  of  the  model 
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which  indicated  a sufficient  margin  of  safety  for  wind  tunnel  testing  at 
this  load  condition. 

The  NASTRAN  vibration  analysis  was  accomplished  with  the  model 
shown  in  Figure  17.  The  model  included  the  internally  calculated  plate 
mass,  inertia  and  stiffness,  the  measured  mass  and  inertia  characteristic 
of  the  sleeve,  and  the  additional  stiffness  of  the  sleeve.  The  mass  and 
inertial  characteristics  of  each  sleeve  section,  given  in  Table  3,  were 
simulated  by  pairs  of  equal  masses,  balanced  about  the  section  center  of 
gravity.  The  additional  torsional  stiffness  of  the  sleeve  was  simulated 
by  rigidly  connecting  the  six  bridge  end  points  local  rotational  degrees 
of  freedom  about  the  spanwise  axis  running  parallel  to  the  wing  leading 
edge.  These  two  modifications  to  the  basic  plate  analytical  model 
effectively  accounted  for  the  sleeve's  influence  on  the  model. 


The  ground  vibration  test  revealed  the  analytical  frequencies  and 
mode  shapes  were  not  accurate.  This  problem  was  solved  by  locating  the 
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cantilevered  root  of  the  analytical  model  1.5  inches  inboard  to  compen- 
sate for  the  mount  flexibility.  The  NASTRAN  calculated  natural  fre- 
quencies are  presented  in  Table  15,  Good  correlation  in  frequency  and 
mode  shapes  was  obtained  for  the  four  plate  models.  The  largest  dis- 
crepancy occurred  with  the  aluminum  plate  model  data  where  the  analytical 
fifth  and  sixth  modes  occurred  in  reverse  order.  This  was  not  considered 
a problem  as  only  the  first  three  modes  were  used  in  the  flutter  and 
divergence  calculations.  The  NASTRAN  mode  shapes  for  each  plate  model 
are  presented  in  Appendix  A. 


TABLE  15 

NASTRAN  ANALYSIS  WIND  TUNNEL  MODEL  NATURAL  FREQUENCIES 


Mode 

A1  uminum 

Non rota ted 

Graphite-epoxy 

7.5°  Rotated 

15°  Rotated 

IB 

3.20 

3.51 

3.37 

3.11 

2B 

16.30 

17.75 

17.10 

15.66 

IT 

36.12 

32.31 

31.73 

31.36 

3B 

42.94 

46.51 

44.91 

41.01 

2T 

85.23 

70.28 

69.23 

69.19 

4B 

82.83 

88.93 

86.09 

78.55 

Units:  Hz 


Once  the  vibration  models  gave  satisfactory  results,  the  flutter 
and  divergence  calculations  were  accomplished.  The  frequency  versus 
velocity  and  damping  versus  velocity  plots  obtained  from  NASTRAN  for 
each  of  the  twenty  plate/sweep  angle  combinations  are  contained  in 
Appendix  C.  In  every  case,  the  first  bending  mode  was  the  mode  that 
diverged.  However,  flutter  occurred  before  divergence  In  some  models  at 
the  zero  sweep  angle.  The  flutter  and  divergence  predictions  from 
NASTRAN  are  presented  In  Tables  16  and  17. 


33 
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NASTRAN  ANALYSIS  FLUTTER 

DYNAMIC  PRESSURES 

! Sweep 

Aluminum 

Nonrotated 

Graphite-epoxy 

7.5°  Rotated 

15°  Rotated 

' 0° 

56.9 

34.4 

36.2 

37.4 

-15° 

66.3 

39.5 

39.8 

42.0 

91.2 

53.2 

50.9 

52.5 

-45° 

139.3 

113.6 

125.8 

126.9 

-60° 

94.2 

88.9 

87.2 

76.7 

Units:  psf 


TABLE  17 

NASTRAN  ANALYSIS  DIVERGENCE  DYNAMIC  PRESSURES 


SECTION  V 


WIND  TUNNEL  TESTING  AND  CORRELATION  OF  RESULTS 

1.  WIND  TUNNEL  TESTING  AND  PROJECTION  METHODS 

The  nature  of  divergence  does  not  allow  testing  near  the  divergence 
velocity  because  there  is  usually  little  chance  of  recovery.  Since 
20  wing  configurations  were  to  be  tested,  each  with  the  same  foam 
sleeve,  it  was  necessary  to  use  a testing  technique  that 'limited  the 
possibility  of  damage  to  the  sleeve,  as  well  as  the  plates  and  wind 
tunnel . 

The  subcritical  divergence  testing  began  at  50 % of  the  analyti- 
cally predicted  divergence  velocity.  At  this  velocity,  the  model 
was  positioned  at  the  angle  of  attack  where  no  bending  strain  was 
observed.  The  model  angle  of  attack  was  increased  in  increments  of  one 
degree,  and  bending  and  torsional  strain  readings  were  recorded  at  the 
nominal  angles  of  attack  of  1 , 2 and  3°.  The  model  was  returned  to  the 
initial  angle  of  attack,  the  tunnel  dynamic  pressure  was  increased, 
and  strain  readings  were  again  taken  at  the  three  angles  of 
attack. 

A minimum  of  six  dynamic  pressures  and  therefore  six  sets  of  strain 
readings  were  recorded  up  to  80%  of  the  projected  divergence  velocity. 
The  divergence  velocity  was  projected  by  two  techniques,  the  divergence 
index  method  and  a Southwell -type  method,  which  were  programmed  in  a 
hand-held  calculator.  As  each  set  of  data  was  read,  it  was  stored  in 
the  calculator  and  the  divergence  velocity  projections  were  updated. 

The  divergence  index  method  was  developed  at  the  NASA  Langley 
Research  Center.  The  discussion  given  here  was  obtained  by  the  authors 
from  Mr.  W.  H.  Reed,  Chief  of  the  Aeroelasticity  Branch  at  the  NASA 
Langley  Research  Center.  For  the  wing  shown  in  Figure  18,  the  lift  of 
the  wing  is  the  sum  of  the  rigid  lift  due  to  angle  of  attack,  a,  and  the 
incremental  lift  due  to  the  angle  of  attack  induced  by  the  fl exibil ity* 

0,  of  the  wing 

L = qSCL  (a  + 0)  (5) 

a 

The  restoring  force  of  the  wing 

Fr  = k0  (6) 
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Figure  18.  Wing  Angle  of  Attack  Components. 

(where  k represents  the  stiffness  of  the  wing  structure)  is  equal  to  the 
wing  lift. 

Equating  Equations  5 and  6,  and  rearranging  terms 


6 = 


qSC^  a 
a 

k-qSC, 


(7) 


"a 


At  divergence  the  denominator  of  Equation  7 is  equal  to  zero.  Thus  the 
restoring  force  is  just  able  to  counter  the  flexible  lift,  giving 

(8) 


a 

and  the  divergence  dynamic  pressure  Is 

n - k 

qD  ' SC[ 


(9) 


a 


Substituting  Equation  8 into  Equation  7 and  solving  for  the  angle  of 
attack  as  a function  of  the  angle  of  attack  Induced  by  the  wing  flexibil- 
ity results  in 

On 

(10) 


Measuring  strain  in  the  wind  tunnel  model  with  a strain  gage  located 
on  the  plate  near  the  root  provides  an  indication  of  the  wing  deflection 
due  to  the  flexible  lift.  The  relation  of  strain  to  deflection  Is 

0 = Be  (11) 

where  B is  a constant.  Substituting  Equation  11  into  Equation  10  results 

in  % 

a = (-~  - 1 )Be  (12) 


‘ f 
i 


T l 
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which  illustrates  that  the  strain  varies  linearly  with  angle  of  attack 
fcr  a constant  dynamic  pressure. 

The  slope  of  the  angle  of  attack  versus  strain  curve  for  the  nth 
dynamic  pressure  is 


A 


n 


de  \ 


03) 


Dividing  the  slope  of  the  first  dynamic  pressure  line  by  the  slope  of 


the  nth  dynamic  pressure  line  results  in  the  equation 

h.  B(vqi : 1} . yqi  - qn/qD 

A„  " B(qr,/q„  - T)  " 1 - 


n " ' "v  ’’n 
The  divergence  index  is  defined  as 

and  when  substituted  into  Equation  14,  results  in 

A _i  % 

n qD 


04) 


(15) 


06) 


Thus,  the  divergence  index  varies  linearly  with  the  dynamic  pressure. 

As  the  dynamic  pressure  approaches  the  divergence  dynamic  pressure,  the 
divergence  index  approaches  zero.  An  example  of  the  use  of  the  diver- 
gence index  is  outlined  in  the  following  paragraphs. 

The  angle  of  attack  versus  strain  data  for  the  nonrotated  model 
at  -15°  sweep  Is  presented  In  Table  18  and  shown  plotted  in  Figure  19, 

The  slope  for  the  first  set  of  data  is  calculated  and  used  as  the  refer- 
ence. The  slope  for  each  subsequent  set  of  data  is  calculated  and  substi- 
tuted into  Equation  15  to  calculate  the  divergence  indices  which  are 
tabulated  in  Table  19. 

The  divergence  Indices  are  plotted  versus  dynamic  pressure  In 
Figure  20.  For  a dynamic  pressure  of  zero  the  divergence  index  Is  one, 
therefore,  the  linear  relation  between  the  divergence  index  and  the 
dynamic  pressure  is  fitted  by  a least  squares  method  and  forced  through 
An  = 1 . The  divergence  dynamic  pressure  is  the  intersection  of  this 
line  and  the  dynamic  pressure  axis.  For  this  example,  the  divergence 
dynamic  pressure  is  projected  to  be  11,9  psf. 
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TABLE  19 

DIVERGENCE  INDEX  PROJECTIONS 
FOR  THE  NONROTATED  MODEL,  A = -15° 


Dynamic 

Pressure 

(psf) 

Divergence 

Index 

A 

Divergence 

Projection 

(psf) 

4.03 

- 

- 

4.73 

.688 

15.1 

5.48 

.520 

12.9 

5.88 

.490 

12.4 

6.29 

.428 

11.9 

6.72 

.432 

11.9 

l 


) 
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The  second  subcritical  projection  method  used  in  the  wind  tunnel 
testing  is  an  adaptation  of  Southwell's  technique  [17].  Southwell's 
technique  was  originally  developed  to  project  beam  buckling  by  reducing 
the  influence  of  geometric  imperfections.  The  similarity  between  wing 
divergence  and  conventional  buckling  of  structures  has  been  noted  by 
many  (Reference  18  is  one  example).  In  discussions  concerning  wind 
tunnel  testing  for  divergence,  it  was  proposed  that  a Southwell  type 
technique  could  be  used  as  a subcritical  projection  method. 

Equation  12  can  be  rearranged  as  follows 

f ■ (0  + B.)  (17) 

For  constant  angle  of  attack,  Equation  17  is  a linear  relationship 
between  e/q  and  e where  the  slope  is  the  inverse  of  the  divergence  dynamic 
pressure: 


^ - Vq, 


Equation  18  is  analogous  to  the  relation  between  load  and  beam  deflection 
in  Reference  17. 

As  in  the  case  of  the  divergence  index  method,  strain  is  measured 
at  each  dynamic  pressure  and  angle  of  attack.  For  two  or  more  dynamic 
pressures,  the  strain  data  at  constant  angle  of  attack  is  fitted  by  a 
least  squares  method.  The  inverse  of  the  slope  of  this  fit  is  the 
projected  divergence  dynamic  pressure.  As  new  strain  data  is  collected, 
the  divergence  dynamic  pressure  is  updated,  as  shown  in  Table  20  for  the 
data  presented  in  Table  18.  The  data  obtained  for  the  nonrotated  model 
is  presented  in  Figure  21.  The  subcritical  projections  compare  favorably 
with  the  divergence  index  projections. 

Each  of  the  20  wing  configurations  was  tested  subcritically 
using  both  the  divergence  index  and  Southwell  methods.  Only  the  projec- 
tions are  presented  in  this  report  since  the  compilation  of  the  interme- 
diate data  is  voluminous.  However,  one  set  of  data,  the  data  for  a 
configuration  that  is  divergence  free,  is  interesting  enough  to  Include 
in  this  report. 

Table  21  is  a tabulation  of  the  dynamic  pressure,  angle  of  attack 
and  strain  data  obtained  while  testing  the  15°  rotated  model  at  0° 
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sweep.  The  slopes  from  the  angle  of  attack  versus  strain  data,  Figure 
22,  were  used  in  the  calculation  of  the  divergence  index.  Figure  23  is 
the  divergence  index  plot  and,  although  there  is  a lot  of  scatter,  the 
trend  is  obvious.  The  divergence  index  for  a divergence-free  wing  is 
greater  than  1.0 

The  Southwell  plot  is  shown  in  Figure  24.  Again,  there  is  scatter 
in  the  slopes,  yet  a trend  exists.  For  a divergence-free  wing  the 
Southwell  plot  slope  is  negative. 

After  subcritical  testing  was  completed  on  the  20  configurations, 
the  aluminum  model  at  -30°  sweep  was  selected  to  obtain  a "hard",  or 
actual,  divergence  data  point.  The  Southwell -type  subcritical  divergence 
projection  technique  was  used  up  to  80%  of  the  divergence  velocity.  The 
wing  angle  of  attack  was  adjusted  to  minimize  the  bending  strain  and 
fixed  for  the  remainder  of  the  run.  The  velocity  of  the  tunnel  was 
raised  incrementally  until  the  wing  divergence  occurred  at  which  time 
the  wind  tunnel  was  immediately  shut  down. 

Fortunately,  divergence  did  not  cause  damage  to  the  model.  As 
previously  described,  the  airfoil  sleeve  was  sectioned  in  order  to 
reduce  the  stiffness.  However,  under  the  large  deflections  associated 
with  the  diverging  model,  the  sleeve  sections  pressed  against  one  another 
causing  the  sleeve  to  restrain  the  model.  This  characteristic  of  the 
model  made  it  possible  to  find  the  actual  divergence  points  for  each  of 
the  four  models.  Several  repeated  cases  demonstrated  that  the  models 
gave  consistent  divergence  results.  The  models  were  not  tested  to 
divergence  at  the  -60°  sweep  position  because  of  the  greater  possibility 
of  damage  due  to  the  sleeve  striking  the  fairing.  Similarly,  no  diver- 
gence points  were  obtained  at  zero  sweep  because  of  the  possibility  of 
encountering  a high  frequency  (~25  Hz)  flutter  instability. 

2.  DISCUSSION  OF  RESULTS  AND  CORRELATION 

The  divergence  dynamic  pressures  measured  during  the  wind  tunnel 
tests  have  been  nondimensionalized  by  plate  weight  divided  by  plate  area 
and  plotted  versus  wing  leading  edge  sweep  in  Figure  25.  For  those 
configurations  where  the  actual  divergence  points  were  not  obtained  the 
Southwell  predictions  were  plotted.  All  of  the  divergence  dynamic 
pressures,  analytical  and  experimental,  are  presented  iri  Table  22. 
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Figure  23.  Divergence  Index  Projection  for  the 
15°  Rotated  Model,  A = 0°. 
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Figure  24.  Southwell  Plot  for  the  15°  Rotated  Model,  A = 0°. 
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Figure  25.  Nondimensional  Divergence  Dynamic  Pressures  versus  Sweep. 

Two  observations  concerning  the  use  of  a composite  material  in 
place  of  aluminum  can  be  made  from  the  comparison  in  Figure  25.  One  is 
that  for  all  sweeps  the  composite  plates  are  more  effective  per  unit 
weight  than  the  aluminum  plate  in  preventing  divergence.  The  second 
observation  is  that  the  divergence  speed  of  the  model  could  be  altered 
by  simply  rotating  the  composite  laminate  in  relation  to  the  reference 
line  of  the  wing.  This  is  especially  evident  at  sweep  angles  between 
0°  and  -20°  where  forward  rotation  of  the  composite  laminate  has  the 
greatest  effect. 

Rotating  the  laminate  is  a form  of  aeroelastic  tailoring.  The 
effect  of  laminate  rotation  may  be  further  appreciated  by  studying 
Figure  26  which  presents  the  laminate  bending  and  torsional  stiffness 
and  the  coupling  parameter  for  each  of  the  composite  plates  as  defined 
by  the  method  of  Reference  5.  Torsional  stiffness  is  nearly  constant 
between  plus  and  minus  5°  rotation  and  increases  sharply  at  higher 
rotation  angles.  At  15°  rotation,  torsional  stiffness  is  about 
60%  higher  than  at  0°  rotation.  The  coupling  parameter  has  a nearly 
constant  slope,  increasing  negatively  from  0°  r-  .ation.  Negative  cou- 
pling parameter  produces  a wash-out,  bend-twist  characteristic  about  the 
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Figure  26.  Stiffness  Variation  Due  to  Rotation  of  [0- , ( -45 ,+45)2]s 
Graphite-epoxy  Laminates. 

! 

structural  axis.  At  the  low  forward  sweeps,  the  divergence  mode  is  j 

primarily  a torsion  mode.  Therefore,  the  increasing  torsional  stiffness  j 

and  decreasing  coupling  parameter  due  to  15°  rotation  have  the  greatest  « 

effect  countering. the  wash-in  tendencies  and  increasing  the  divergence  | 

speed.  | 

Figure  25  shows  that  the  15°  rotated  model  has  the  lowest  diver- 
gence dynamic  pressure  of  the  three  composite  models  at  -60°  sweep. 

At  -45°  sweep,  the  15°  rotated  model  has  a divergence  dynamic  pressure 

that  is  less  than  the  7.5°  rotated  model.  At  the  greater  forward  sweeps, 

the  divergence  mode  is  primarily  bending  and  bending  stiffness  becomes 

predominant  in  determining  divergence  dynamic  pressure.  Figure  26  shows  i 

the  bending  stiffness  is  nearly  constant  over  the  ±5°  rotation  range, 

but  is  about  11%  less  for  15°  rotation  than  for  the  nonrotated  laminate. 

Although  the  torsional  stiffness  and  wash-out  coupling  are  greatest  for 
this  laminate,  the  bending  stiffness  is  lowest  resulting  in  low  diver- 
gence speeds  at  the  greater  forward  sweeps. 
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Figure  27  presents  the  "hard"  divergence  points  compared  with  the 
subcritical  projections  based  on  the  Southwell  method  for  the  7.5° 
f rotated  model.  Similar  results  were  obtained  for  the  c.her  models.  The 

Southwell  method  projected  divergence  dynamic  pressures  within  10°/  for 
the  configurations  where  subcritical  data  were  obtained  at  test  points 
j greater  than  50%  of  the  divergence  dynamic  pressure.  It  was  not  possible 

I l to  obtain  data  at  50%  of  the  divergence  dynamic  pressures  for  the  0° 

f sweep  cases  due  to  the  low  flutter  speeds.  Therefore,  the  quality  of 

convergence  of  the  projections  at  0°  sweep  was  poor. 
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Figure  27.  Comparison  of  Measured  and  Projected  Divergence 
Dynamic  Pressures  for  the  7.5°  Rotated  Model. 


As  seen  In  Figure  28  for  the  aluminum  model,  the  divergence  dynamic 
pressure  decreases  rapidly  when  the  wing  Is  swept  from  0°  to  -15°  and 
remains  nearly  constant  from  -30°  to  -60°  sweep.  This  trend  Is  predicted 
very  well  by  TS0  and  NASTRAN  analyses.  The  TS0  and  NASTRAN  analytical 
predictions  using  doublet  lattice  aerodynamics  are  In  close  agreement 
with  the  experimental  data.  As  discussed  in  a previous  section,  the  TS0 
analysis  with  Woodward  aerodynamics  consistently  predicts  slightly 
higher  divergence  dynamic  pressures  than  the  TS0  analysis  with  doublet 


lattice  aerodynamics. 


Figure  28.  Comparison  of  Analytical  and  Test  Divergence  Dynamic 
Pressures  for  the  Aluminum  Model. 


The  analytical  and  test  divergence  data  for  the  nonrotated  model 
: are  presented  in  Figure  29.  For  the  0°  sweep  case,  the  flutter  dynamic 

i pressure  was  lower  than  the  divergence  dynamic  pressure.  Consequently, 

testing  was  restricted  to  below  35%  of  the  divergence  dynamic 
pressure,  and  the  subcritical  projections  did  not  converge.  The  corre- 
,■  lation  between  test  and  analytical  divergence  dynamic  pressures  at  -15° 

sweep  Is  poor.  This  poor  correlation  caused  concern,  and  thus,  vibra- 
tion and  load-deflection  tests  were  performed  after  the  wind  tunnel 
tests.  The  results  and  a discussion  of  the  results  are  presented  later 
in  this  section  of  the  report. 

Figure  30  presents  the  analytical  and  test  divergence  data  for 
the  7.5°  rotated  model.  Rotation  of  the  laminate  7.5“  forward  of  the 
reference  line  significantly  Increases  the  divergence  dynamic  pressure 
l at  0°  sweep.  As  was  the  case  for  the  nonrotated  composite  plate  at  this 

j,  sweep,  the  Southwell  divergence  projection  did  not  converge  because 

| it  was  not  possible  to  test  to  sufficiently  high  dynamic  pressures  due 

\ to  the  low  flutter  speeds.  The  calculated  divergence  dynamic  pressure 

j at  0°  sweep  for  the  7.5°  rotated  model  was  at  least  four  times  the 
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LEADING  EDGE  SWEEP 

Figure  29.  Comparison  of  Analytical  and  Test  Divergence  Dynamic 
Pressures  for  the  Nonrotated  Model . 


LEADING  EDGE  SWEEP 

Figure  30.  Comparison  of  Analytical  and  Test  Divergence  Dynamic 
Pressures  for  the  7.5°  Rotated  Model. 
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divergence  dynamic  pressure  for  the  nonrotated  model  (Table  22).  The 
NASTRAN  and  TSO  doublet  lattice  analytical  results  compare  very  well  with 
test  data  for  this  plate  at  sweeps  greater  than  15°  forward.  The  TSO- 
Woodward  analytical  results  are  consistently  higher  than  the  test  data. 

Figure  31  presents  the  analytical  and  test  divergence  data  for  the  • 
15°  rotated  model.  As  seen  in  this  figure  and  in  Table  20,  the  analyti- 
cal divergence  dynamic  pressure  is  very  high  at  0°  sweep.  As  discussed 
earlier,  the  subcritical  techniques  indicated  a divergence-free  wing  for 
this  plate  and  sweep.  At  -15°  sweep,  divergence  was  very  difficult,  if 
not  impossible,  to  define.  Subcritical  projections  were  obtained,  but 
as  the  tunnel  dynamic  pressure  was  increased,  the  projected  divergence 
dynamic  pressure  increased.  An  unusual  phenomenon,  a low  amplitude,  low 
frequency  (0.5  Hz)  sinusoidal  oscillation  of  the  wing,  occurred  at 
approximately  the  analytically  predicted  divergence  speed.  As  the 
dynamic  pressure  increased,  the  amplitude  of  the  oscillation  increased 
while  the  frequency  remained  constant.  The  cause  of  the  phenomenon  and 
its  mechanism  are  not  understood.  It  is  not  predicted  by  any  of  the 
analytical  methods  and  thus  may  be  associated  with  a characteristic  of 
the  model  tested.  Thus  a meaningful  comparison  between  the  analytical 
and  test  results  is  not  possible.  At  -30°  and  -45°  sweep,  comparisons 
of  analytical  and  test  results  are  good.  However,  at  -60°  sweep, 
the  Southwell  projection  is  12.4%  less  than  the  lowest  analytical 
prediction.  Also,  the  trend  of  increasing  divergence  dynamic  pressure 
predicted  by  the  analysis  for  sweeps  greater  than  45°  forward  is  opposite 
to  the  trend  observed  in  the  test  results. 

The  divergence  characteristics  varied  with  sweep  angle  and  structur- 
al plate.  Generally,  the  severity  of  divergence,  described  as  the  rate 
of. change  of  deflection  as  the  wing  diverged,  was  greater  at  the  higher 
forward  sweep  angles.  Rotating  the  composite  laminate  forward  lessened 
the  severity  of  the  divergence  at  all  forward  sweep  angles  where  "hard" 
divergence  points  were  obtained.  At  -15°  sweep,  the  rate  of  deformation 
associated  with  divergence  was  mild,  while  at  -45°  sweep,  the  rate  was 
rapid.  The  rapid  rate  of  deformation  caused  the  sleeve  sections  to  com- 
press as  a spring  which  resulted  in  a post-divergence  oscillation  (7  Hz). 

At  -60®  sweep,  the  test  results  are  consistently  lower  than  all  of 
the  analytical  results  for  all  models.  Previous  testing  involving  the 
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LEADING  EDGE  SWEEP 

Figure  31.  Comparison  of  Analytical  and  Test  Divergence  Dynamic 
Pressures  for  the  15°  Rotated  Model. 

fairing  revealed  that  turbulence  was  generated  around  the  cavity  of  the 
fairing.  For  this  test,  cover  plates  were  used  to  minimize  the  cavity 
and  reduce  the  turbulence.  The  size  of  the  cavity  was  largest  at  -60°. 

The  turbulence  generated  by  the  fairing  cavity  could  have  affected  the 
aerodynamic  loading  on  the  inboard  aft  portion  of  the  model.  Thus  the 
center  of  pressure  would  be  more  forward.  Since  the  analysis  does  not 
account  for  the  turbulence  near  the  fairing,  the  calculated  center  of 
pressure  would  be  further  aft  than  the  actual  location  on  the  model. 

For  -60°  sweep,  the  effect  of  cavity  turbulence  would  be  greatest, 
possibly  explaining  the  difference  between  the  analytical  and  test  results. 

The  effect  of  laminate  rotation  on  loading  is  illustrated  In  Figure 
32.  Measured  strain  is  plotted  versus  dynamic  pressure  for  the  nonrota- 
ted  and  15°  rotated  models  for  -30°  sweep  at  3°  angle  of  attack. 

For  dynamic  pressures  greater  than  50%  of  the  divergence  dynamic 
pressure  of  the  nonrotated  model,  the  strain  level  is  lower  for  the  15° 
rotated  model  than  for  the  nonrotated  model.  Therefore,  increasing  the 
divergence  dynamic  pressure  by  laminate  rotation  decreased  the  level  of 
strain  under  aerodynamic  loading. 
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Figure  32.  Comparison  of  Strain  Levels  for  Two  Composite  Models. 

3.  POST  TUNNEL  TESTING  AND  ANALYSIS 

As  indicated  in  the  previous  section,  a comparison  of  the  test 
results  with  analytical  results  shows  some  inconsistencies.  Eliminating 
the  unswept  results  from  this  discussion  due  to  the  unreliable  projec- 
tions because  of  flutter,  and  the  -60°  sweep  positions  due  to  the  ques- 
tion  of  the  analytical  aerodynamics  at  this  high  sweep,  the  following 
discussion  will  concentrate  on  the  three  interim  sweep  positions  (-15°, 
-30°  and  -45°). 

Comparing  the  analytical  results  with  the  hard  divergence  points, 
the  aluminum  and  the  7.5°  rotated  models  differ  by  less  than  5%.  How- 
ever, analysis  of  the  nonrotated  model  predicted  unconservative  diver- 
gence dynamic  pressures;  in  the  case  of  the  15®  forward  sweep,  the 
analysis  predicted  a 24%  higher  divergence  dynamic  pressure  than 
was  measured.  For  the  15®  rotated  model,  the  analytical  results  are 
conservative;  divergence  at  15®  forward  sweep  was  undefinable. 


During  manufacturing  of  the  sleeve,  the  channels  were  mounted  on 
the  aluminum  plate  which  was  used  as  part  of  the  mold.  When  the  sleeve 
was  attached  to  the  aluminum  plate,  the  model  was  in  an  unstressed 
state.  However,  the  composite  plates  are  thinner,  requiring  the  sleeve 
channels  and  dowels  to  compress  the  foam  when  bolted  to  the  plate.  The 
resulting  tension  in  the  bolts  caused  stress  to  be  applied  to  the  com- 
posite plates.  This  stress  has  the  general  effect  of  stiffening  the 
model,  particularly  in  torsion,  and  increasing  the  torsional  frequencies. 

Prior  to  the  wind  tunnel  testing,  the  sleeve  never  experienced  a 
high  load  condition.  During  wind  tunnel  testing,  some  model  configura- 
tions were  loaded  so  the  deflection  was  approximately  a foot  at  the  tip. 
While  testing  the  aluminum  model  to  find  the  "hard"  divergence  points, 
the  model  experienced  higher  deflections  and,  in  some  cases,  severe  high 
amplitude  oscillations.  The  result  of  the  high  deflection  was  to  loosen 
the  bond  in  the  sleeve  between  the  channels  and  the  foam.  In  order  to 
provide  a better  correlation  between  the  analyses  and  test,  the  ground 
vibration  tests  were  repeated  and  the  analytical  models  were  redefined. 

Table  23  compares  the  results  from  the  ground  vibration  tests.  The 
greatest  change  is  the  reduction  in  the  first  torsional  frequency. 

Since  the  mass  of  the  model  remains  unchanged,  the  model  torsional 
stiffness  must  have  been  lowered.  Although  not  shown,  the  plate  fre- 
quencies were  unchanged.  Therefore,  the  reduction  in  model  torsional 
stiffness  is  due  to  the  reduction  in  stiffness  due  to  the  sleeve. 

Input  to  NASTRAN  and  TSO  analyses  were  modified  to  account  for 
changes  in  the  composite  models.  The  stiffness  of  the  finite  element 
beams  in  TSO  and  the  constrained  plate  elements  in  NASTRAN  was  reduced 
until  the  analytical  model  frequencies  matched  the  frequencies  measured 
during  the  ground  vibration  test  conducted  after  the  wind  tunnel  test. 

The  results  of  the  revised  analyses  are  shown  in  Table  22.  As  expected, 
the  nonrotated  post  test  analyses  showed  an  overall  reduction  in  the 
divergence  speeds.  The  revised  analyses  predict  divergence  dynamic 
pressures  less  than  8%  above  the  test  values.  For  the  7.5°  rotated 
model,  the  analytical  results  changed  very  little  (<1%).  For  the  15° 
rotated  model,  the  revised  analyses  predicted  divergence  dynamic  pres- 
sures that  are  within  4%  of  the  test  results. 
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TABLE  23 

COMPARISON  OF  MODEL  NATURAL  FREQUENCIES 
MEASURED  BEFORE  AND  AFTER  THE  WIND  TUNNEL  TEST 


Graphite-epoxy 


Mode 

Aluminum 

Non rota ted 

7.5°  Rotated 

15°  Rotated 

0) 

(2) 

(1) 

(2) 

(1) 

(2) 

(1) 

(2) 

X 

IB 

3.32 

3.29 

3.57 

3.51 

3.40 

3.39 

3.19 

3.18 

% 

2B 

16.97 

16.78 

17.04 

17.67 

16.96 

16.93 

15.36 

15.40 

IT 

36.09 

36.95 

31.57 

26.90 

32.04 

26.58 

31 ,42 

27.76 

3B 

43.48 

43.61 

46.53 

45.68 

44.36 

43.48 

40.70 

40.20 

y 

2T 

87.07 

88.87 

64.88 

63.16 

65.88 

65.85 

68.74 

71.68 

} 

1 

4B 

81  .85 

81.43 

86.22 

84.28 

83.08 

81.91 

75.99 

76.28 

i 

5 

Units:  Hz 

GVT  performed  before  wind  tunnel  test. 
GVT  performed  after  wind  tunnel  test. 


Although  the  overall  comparison  has  been  improved  by  matching  the 
analysis  to  the  post  wind  tunnel  testing,  the  increase  of  divergence 
dynamic  pressure  for  the  15°  rotated  model  was  not  expected  since  the 
model  torsional  stiffness  had  decreased.  A possible  explanation  for 
the  increase  in  divergence  dynamic  pressure  is,  as  the  plies  are  rotated 
forward,  the  coupling  between  the  sleeve  and  plate  increases.  For  the 
15°  rotated  model,  the  coupling  caused  by  the  sleeve  must  have  been 
detrimental,  that  is,  the  coupling  caused  a reduction  in  the  divergence 
dynamic  pressure.  Therefore,  a reduction  in  the  sleeve  stiffness  would 
cause  an  increase  in  the  divergence  dynamic  pressures.  For  the  nonro- 
tated  model,  the  sleeve  channels  are  nearly  perpendicular  to  the  primary 
bending  plies  and  provide  minimal  coupling  between  bending  and  torsion. 
Therefore,  a reduction  in  the  torsional  stiffness  component  would  result 
in  a reduction  of  the  divergence  dynamic  pressure.  This  is  especially 
evident  at  low  forward  sweeps  where  torsional  stiffness  has  the  greatest 
effect. 
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SECTION  VI 
CONCLUSIONS 


I 

\ 


Results  of  analyses,  laboratory  tests,  and  wind  tunnel  tests  of  a 
rather  simple,  variable  sweep  model  that  could  be  swept  forward  15°, 

30°,  45°  and  60°  from  the  leading  edge  unswept  position,  and  could  also 
accommodate  structural  plates  of  aluminum  and  graphite-epoxy  materials, 

* have  illustrated  the  principle  involved  In  the  structural  design  tech- 

nology of  aeroelastic  tailoring.  Several  conclusions  can  be  drawn  from 
this  research  which  are  briefly  discussed  below. 

The  high  stiffness  to  weight  ratio  of  the  graphite-epoxy  is  more 
efficient  than  aluminum  in  providing  the  stiffness  required  to  increase 
the  divergence  dynamic  pressure.  The  ability  to  tailor  the  composite 
material  for  bend-twist  coupling  significantly  adds  to  the  efficiency 
from  a weight  standpoint. 

It  has  been  shown  that  simply  rotating  a 0°±45°  composite  laminate 
forward  significantly  increases  the  divergence  dynamic  pressure  of  a 
forward  swept  wing  at  leading  edge  sweep  angles  to  about  -20°.  A 
smaller  increase  in  divergence  dynamic  pressure  occurs  at  -30°  and  -45° 
sweep.  A reversal  in  the  trend  occurs  at  -60°  sweep. 

The  analytical  methods  used  predict  the  divergence  dynamic  pressures 
very  well  for  all  models  at  -30°  and  -45°  sweep.  For  -15°  sweep,  the 
correlation  between  analytical  and  test  results  is  very  good  for  all 
models  except  the  15°  rotated  model  where  an  unusual  oscillatory  phenom- 
enon occurred  in  the  wind  tunnel.  The  large  deflections  obtained  during 
testing  apparently  caused  a reduction  in  torsional  stiffness  of  the 
composite  models,  probably  due  to  loosening  of  the  aerodynamic  sleeve. 
This  effect  was  most  noticeable  at  the  0°  sweep  angle  and  affected  the 
correlation  of  analytical  and  test  results.  The  analytical  results 
were  consistently  higher  than  the  measured  values  at  -60c  sweep,  appar- 
ently associated  with  turbulence  from  the  fairing  cavity. 

The  divergence  dynamic  pressures  predicted  by  the  static  and 
dynamic  (velocity-damping)  analyses,  using  Woodward  and  Doublet  Lattice 
aerodynamics,  respectively,  agree  favorably  at  sweep  angle  of  -15°,  -30° 
and  -45°.  The  Woodward  static  analysis  is  least  accurate  at  -60°  sweep, 

* predicting  higher  dynamic  pressure  than  measured. 
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While  the  divergence  dynamic  pressure  is  not  significantly  increased 
at  -30°  sweep  by  rotating  the  composite  laminate  15°  forward,  the  load 
level  at  a given  angle  of  attack  is  significantly  reduced  at  dynamic 
pressures  greater  than  60%  of  the  divergence  dynamic  pressure. 

The  subcritical  projection  methods  described  herein  accurately 
predict  divergence  dynamic  pressure  at  80%  and  less  of  the  divergence 
dynamic  pressure.  It  may  be  possible  to  use  methods  like  these  in 
flight  testing  for  divergence. 

The  severity  of  the  motion  of  the  wing  at  divergence  onset  Increases 
with  forward  sweep. 


APPENDIX  A 


MODE  SHAPES 

Three  sets  of  mode  shapes  are  presented  in  the  following  pages. 
Figures  A-l  through  A-4  are  the  experimentally  measured  mode  shapes. 
Figures  A-5  through  A-8  are  the  mode  shapes  resulting  from  the  ISO 
analysis  of  each  model.  Figures  A-9  through  A-l 2 are  the  mode  shapes 
resulting  from  the  NASTRAN  analysis  of  each  model. 
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b.  Second  Bending 
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b.  Second  Bending 


Figure  A-8.  TSO  analysis  mode  shapes  for  the  15°  rotated  model. 


Figure  A-9.  NASTRAN  analysis  mode  shapes  for  the  aluminum  model. 
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Figure  A-10.  NASTRAN  analysis  mode  shapes  for  the  nonrotated  model. 
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Figure  A-ll.  NASTRAN  analysis  mode  shapes  for  the  7.5°  rotated  model 


APPENDIX  B 


STACKING  SEQUENCE  ANALYSIS 

The  tables  presented  in  this  Appendix  represent  some  of  the  quali 
tative  analyses  performed  with  the  CWING  computer  procedure.  The  aero 
dynamics  were  calculated  for  the  airfoil  planform.  The  divergence 
dynamic  pressures  were  calculated  for  the  aluminum  plate  and  eight 
graphite-epoxy  plates.  Each  graphite-epoxy  plate  had  plies  of  0°  and 
±45°  in  varying  percentages  ranging  from  all  0°  to  all  ±45? 
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APPENDIX  C 
V-g  AND  V-w  CURVES 


Two  sets  of  V-g  and  V-w  curves  are  presented  in  the  following 
pages.  Each  set  is  comprised  of  one  V-g  and  one  V-w  for  each  of  the 
20  wing  configurations.  The  first  set  resulted  from  the  TSO  analysis 
and  the  second  set  resulted  from  the  NASTRAN  analysis.  All  the  results 
presented  in  this  Appendix  were  calculated  for  sea  level  density. 
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Figure  C-2.  TSO  analysis  V-g  and  V- 
Aluminum  model , A = -15 


°e 


TSO  analysis  V-g  and  V-to  curves. 
Aluminum  model,  A ■ -45°. 


Figure  C-6. 
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Figure  C-ll.  TSO  analysis  V-g  and  V-u>  curves 
7.5°  rotated  model,  A = 0°, 
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Figure  C-13.  TSO  analysis  V-g  and  V-u>  curves 
7,5°  rotated  model,  A = -30°. 
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Figure  C-16 
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Figure  C-17.  TSO  analysis  V-g  and  V-w  curves 
15°  rotated  model,  A = -15°. 
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Figure  C- 18.  TSO  analysis  V-g  and  V-w  curves- 
15*  rotated  model,  A = -30°. 
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Fiqure  C-21 . NASTRAN  analysis  V-g  and  V-w  curves 
Aluminum  model,  A - 0 . 
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Figure  C-22.  NASTRAN  analysis  V-g  and  V-w  curves 
Aluminum  model,  A = -15°. 
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Figure  C-24.  NASTRAN  analysis  V-g  and  V-w  curves 
Aluminum  model,  A = -45°. 
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Figure  C-26.  NASTRAN  analysis  V-g  and  V-ui  curves. 
Nonrotated  model,  A = 0°. 
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Figure  C-32.  NASTRAN  analysis  V-g  and  V-u>  curves 
7.5°  rotated  model,  A = “15°. 
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Figure  C-33.  NASTRAN  analysis  V-g  and  V-w  curves 
7.5”  rotated  model,  A = -30°. 
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Figure  C-3B,  NASTRAN  analysis  V-g  and  V- 1 
15*  rotated  model,  A « -30® 
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APPENDIX  D 
NASTRAN  INPUT  DATA 


The  data  presented  in  this  Appendix  is  the  list  of  input  required 
to  perforin  a dynamic  analysis  using  NASTRAN  of  the  nonrotated  model  at 
-30°  sweep. 
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98-  FLFACT  3 .75  .5  .45  .4  .35  .30  .25  +FL3 

99-  +FL3  0.20  0.15  0.10  0.075  0.050  0.025  0.01  0.0 

100-  FLUTTER  30  K 1 2 3 S 3 

101-  GRDSET  1 6 


in  in  m m in  in  in  unn  ir/  m in  lo  in 

in  in  in  in  in  in  in  in  in  in  m m <\/  (\i  t\j  r>.  r>.  in  in  in  w t\i  i-J  f--.  f>>  in  m u?  im  cm  t\i 

inminMcsjrjNrsi^inininf'Jivicji^r^r-ininin  • • « • * . ♦ • * ♦ • . . . . * * • • • * 

. . . . . . . . . » * . * • * . • v . . ^r-,_,^(v)r0r0u^inii?Lntr)U,r*-rv.r>.aicifn 

• 4 #i-«f-9^n^Y5nrrj<v)MLrttnmr*Hrs»r^r'. 


tq  cocJNin^cn^iiniD^rxtvj^^tntitsiwiDfnMMinooat^niomcotvjcNj^cocof-ininncvj^iDco 
uiineoNi-ta^tnhvPs&icoaJtriDStmfocorj^fOTrvjcM^^oicotfilSjGt^tncvi^iincncoio^r'v^rom*-* 
‘^•-•inMcr»(ncvjfN.oiiMrNW*-<<<4n^(\JwtiM^iT>iqcocncoM-corN’^ff4infeqr>.<nr^vooJioiolS}Mincotq 


fcj^Nn^in»DN®ff'^'-4Nrt,itintONeo^iR^M(*)^u)ioNfl)rt^«(v/<»)Rrin 
rtpjpjRfwiaNcacftrtWHHMM-*  <-t-4<-4NCsiC4M>j(Nirj(M(Mrvip;r)r,x«7nn(vic9o^C7«-^-^^<<‘<>^ 


oooooQOOoaoiDoaooaoiuoooiaooooejociooacicicioooiriaoooci 


ozaz&aiaZat.&.ciZcCQZai.BdQC&.otQ&QZQitxLci'QZat.ai.ataLDt.ccccae.nZciQtQCQ'ctczQZ&.ad.Qt.aiiitaiaiia: 


appendix  e 


MEASURED  INFLUENCE  COEFFICIENT  MATRICES 

The  tables  presented  In  this  Appendix  ere  the  measured  influence 
coefficient  matrices.  The  influence  coefficient  test,  setup  and  data 
reduction  are  presented  in  detail  in  Reference  19. 
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